Hydrodynamic event-plane correlations in Pb+Pb collisions at y / s = 2.76ATeV 
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The recently measured correlations between the flow angles associated with higher harmonics in 
the anisotropic flow generated in relativistic heavy-ion collisions are shown to be of hydrodynamic 
origin. The correlation strength is found to be sensitive to both the initial conditions and the shear 
viscosity of the expanding fireball medium. 

PACS numbers: 25.75.-q, 12.38.Mh, 25.75.Ld, 24.10.Nz 



o- 
in 



(N 
> 
O 
O 
(N 

od 
O 
(N 



1. Introduction. Due to the fluctuating positions of 
the nucleons inside the colliding nuclei at the point of 
impact [l[ and to quantum fluctuations of the quark and 
gluon fields inside those nucleons 0-0] , the initial density 
profiles of the fireballs created in relativistic heavy-ion 
collisions fluctuate in size, shape and magnitude from 
event to event, even for collisions with identical impact 
parameters. One way to characterize these fluctuating 
initial profiles is through a set of harmonic eccentricity 
coefficients e n with associated "participant plane" angles 
$n (see, e.g., [H): 
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where e(r, </>) is the initial energy density distribution in 
the plane transverse to the beam direction at the collision 
point z=0. Heavy-ion collision experiments at the Rela- 
tivistic Heavy Ion Collider (RHIC) and the Large Hadron 
Collider (LHC) have established an extensive data base 
on corresponding anisotropies in the final momentum dis- 
tributions of the emitted charged hadrons, characterized 
by anisotropic flow coefficients v n and their associated 
flow (or "event plane") angles [§-[lo|: 
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Extensive and systematic studies over the last few years 
(for recent reviews see [ll[) have shown that all of the 
measured flow anisotropies can be understood qualita- 
tively, and to a large extent even quantitatively, as hydro- 
dynamic response to the above-mentioned initial-state 
fluctuations, and that the relationship between the har- 
monic flow coefficients v n and corresponding eccentrici- 
ties e n can be used to determine [L2| - [l5| the specific shear 
viscosity (?7/s)qgp of the quark-gluon plasma (QGP) 
that (at RHIC and LHC energies) makes up the fireball 
medium through roughly the first half of its expansion 
history The specific shear viscosity (??/s)qgp turns 
out to be surprisingly small, (?7/s)qgp = (D(l/4n) fl2l - 
[l5j . making the QGP an almost perfect fluid. 

A complete theoretical analysis of all measured har- 
monic flow coefficients v±, . . . , v§ [HI, [l7l - [l9j is not yet 



available since technical shortcuts that for vi and V3 al- 
low one to sidestep the need for event- by- event hydrody- 
namical evolution of large numbers of fluctuating initial 
conditions (lEj fail for the higher-order harmonics [Ioj [. 
making their calculation numerically costly. As pointed 
out in [2(| [U, a complete understanding of the entire 
spectrum of harmonic flow coefficients v n is expected to 
yield strong constraints on the initial conditions and dy- 
namical evolution of heavy-ion collisions, in particular 
the transport coefficients of the fireball medium. The 
authors of 0, [22|-|26| added that correlations between 
the event plane angles fy n of different harmonic order 
can yield valuable additional insights into the initial con- 
ditions. Such correlations were recently measured with 
good precision by the ATLAS Collaboration in Pb+Pb 
collisions at the LHC (27|- We here demonstrate that 
some of the measured final-state event plane correlations 
have a qualitatively different ccntrality dependence from 
the corresponding initial-state participant plane correla- 
tions, and that this characteristic change between initial 
and final state is correctly reproduced by hydrodynamic 
evolution. This provides additional strong support for 
the validity of the hydrodynamic paradigm in relativis- 
tic heavy-ion collisions. Furthermore, we show that the 
measured event-plane correlations are not only sensitive 
to the initial conditions, but also to the shear viscosity of 
the hydrodynamic medium, thus providing an indepen- 
dent constraint for this key transport coefficient. 

2. Methodology. We have used the (2+l)-dimensional 
code VISH2+1 [28| to evolve fluctuating initial energy 
density profiles for Pb+Pb collisions at v / i=2.76^4TeV 
event by event with viscous hydrodynamics. To ex- 
plore the sensitivity to model uncertainties in the ini- 
tial state, we have evolved events from two sets of ini- 
tial conditions obtained from the Monte-Carlo Glauber 
(MC-GhV) and the Monte-Carlo KLN (MC-KLN) mod- 
els [1^, |3(| . We divided each set into centrality classes 
according to the number 7V part of wounded nucleons; for 
each centrality class, we evolved 11,000 events from each 
of the two models. Model parameters were tuned to 
reproduce the pt spectra and elliptic flows of uniden- 
tified charged particles and identified hadrons, as re- 
ported in [lot l3l| . This results in a specific shear vis- 
cosity 77/5=0.08 for MC-Glauber initial conditions and 
the larger value 7//s=0.2 for MC-KLN initial conditions. 
Both the QGP phase and the hadronic phase are evolved 
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FIG. 1: (Color online) Two-plane correlations (cos(jfc($ n — $ m ))), where j is an integer and k is the least common multiple 
(LCM) of n and m [UlU, between pairs of participant-plane angles $„, m for the harmonics (n, m) and multipliers j studied 
in Ref. [23] - Solid (dashed) lines show results for initial density profiles obtained from the MC-Glauber (MC-KLN) model. 




FIG. 2: (Color online) Similar to Fig. [T] but for the corresponding final-state event-plane angles 9„ tm . Filled circles show the 
experimental values measured by ATLAS [27| . The MC-Glauber (solid) and MC-KLN (dashed) initial profiles used in Fig. [1] 
were propagated individually using viscous hydrodynamics with 77/s=0.08 and 0.2, respectively. 



hydrodynamically; particle momentum distributions are 
calculated with the Cooper- Frye prescription, taking into 
account strong decays of all hadron resonances with 
masses up to 2.25 GeV. (We found, however, that the 
event-plane correlations discussed below are almost iden- 
tical for all particle species, so including resonance decays 
is not essential for the present work.) From the resulting 
charged hadron distribution we calculate for each event 
the flow angles according to 

v gm *„ = Jo.5<H<2. 5 d ^ / Pmin PTdPT dcj> p e in ** g^gg d ^ 

/o.6<|„|<2.5 d 1 4™ P T<1 P T d ^P d VP rdp; d <P P 
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employing the same pseudorapidity range 0.5 < \r)\ < 2.5 
and lower px cutoff px > Pmin = 0.5 GeV as used in the 
experimental analysis From these event plane an- 



The ATLAS results were obtained with two independent meth- 
ods: (a) using a calorimetric measurement of transverse energy 
Ex over rapidity range 0.5 < \rj\ < 4.8, and (b) using charged par- 
ticle tracks with px > 0.5 GeV and 0.5 < \rj\ < 2.5. The data from 
method (a) have better precision but are fully compatible with 
those from method (b), within error bars. Since we cannot sim- 
ulate the calorimetric response of ATLAS theoretically, we com- 
pute the event-plane correlations according to method (b), but 
compare them in the figures to the more precise data obtained 
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gles we compute for each event cos(ki^/ ni + . . . +k m ^f nm ) 
for the two-plane (m=2) and three-plane (m=3) correla- 
tions listed in Tables 1 and 2 of Ref. [27| and shown in 
the figures below, and then average this quantity over all 
events in the given centrality class. We compare these 
event-plane correlations with the corresponding correla- 
tions between the initial-state participant plane angles 
calculated from the initial energy density profile of 
each propagated event according to Eq. (TIJ and then av- 
eraged over events in a similar way. 

3. Results. Figures [1] and [2] show the initial and fi- 
nal state two-plane correlations, for the eight different 
combinations of angles and weight factors explored by 
the ATLAS experiment 27| . Each correlation function is 



plotted against collision centrality, with peripheral colli- 
sions (small Ap ar t values) on the left and central collisions 
(large A part ) on the right. Fig. [1] shows that several of 
these correlations are quite sensitive to the model used to 
generate the initial energy density profiles (MC-Glaubcr 
vs. MC-KLN). These model differences in the initial state 
manifest themselves in corresponding model differences 
between the final-state event-plane correlations shown in 
Fig. [2J but they are additionally modified by the different 
shear viscosities r\j s (0.08 and 0.2, respectively) used to 
evolve the initial conditions from the two models. This 
is most clearly seen in the "3-6 correlation", where the 
two models give almost identical initial-state participant- 
plane correlations (cos(6( ( E ) 3— ^e))) (second lower panel 
from the left in Fig. [T]) whereas the corresponding final- 
state event-plane correlators (cos(6(\E , 3— ^e))) exhibit 
significant model differences. This demonstrates the sen- 
sitivity of these event-plane correlations to the specific 
shear viscosity of the expanding fireball medium. 

It is worth emphasizing that several of these two-plane 
correlators exhibit dramatically different centrality de- 
pendences for the initial-state participant-plane and the 
final-state event-plane angles (see, for example, the up- 
per left, two upper right and second lower left panels in 
Figs. Q] and [5]). The difference is largest in peripheral col- 
lisions (small Ap ar t). We believe that this effect is caused 
by a dynamical rotation of the event-plane angles during 
the hydrodynamic evolution [32| , driven by large elliptic 
flow in non-central collisions which leads to mode cou- 
pling between the angles and & n ±2k (where k is an 
integer and the largest coupling coefficient should corre- 
spond to k~l). 2 More detailed analyses will be required 



from method (a). 
2 This is different from the mode-coupling at freeze-out [33ll caused 
by an elliptic (quadrupolc) deformation of the collective flow ve- 
locity appearing in the exponent of the Boltzmann factor in the 
Cooper-Fryc expression for the final particle momentum distri- 
bution that couples v n with v n ±2k- In contrast, in the presence 
of strong elliptic flow the non-linear hydrodynamic evolution be- 
fore freeze- out leads to mode-coupling between the modes n and 
n±2fc (fc integer) for the entire complex flow vector on the left 
hand side of Eq. l(2}. We have checked [32j that the event-plane 
correlations among the finally emitted particles in Figs. [2] and 



to confirm this conjecture. Event-plane rotation in mid- 
central and peripheral collisions has been suggested as 
the mechanism for the decorrelation of the final event- 
plane angle \P n from the initial participant-plane angle 
<&„ observed in Ref. and this interpretation is con- 
sistent with the results of Ref. [34| . Being driven by mode 
coupling, this decorrelation depends strongly on the den- 
sity and shape fluctuations in the initial state (which, 
after all, are the cause for non-zero odd harmonics in 
the final anisotropic flow); we have noticed [32[ that the 
decorrelation is weaker for a viscous fluid than in ideal 
fluid dynamics, reflecting the dynamical smoothing by 
shear viscosity of initial state density fluctuations and 
the flows generated by them, but unambiguous evidence 
for a dynamical rotation of the event planes as the source 
of this decorrelation still needs to be established. 

Figures [3] and [4] show a number of three-plane cor- 
relations studied by the ATLAS experiment [27| . with 
the initial-state participant-plane correlators plotted in 
Fig. [3] and the corresponding final-state event-plane cor- 
relators in Fig. [4j together with the experimental data. 
Again, we observe characteristic sign changes between 
several of the initial-state correlations and their corre- 
sponding final-state correlators. Even if neither of the 
two initial-state models (MC-Glauber and MC-KLN) re- 
produces the experimental data exactly, we find it im- 
pressive that the hydrodynamic model reproduces all the 
qualitative features of the centrality dependences of the 
14 different measured event-plane correlation functions 
correctly: where the data show strong (weak) correla- 
tions, the same is true for the theoretical results, and 
where the data show correlations that increase (decrease) 
from peripheral to central collisions, the same holds for 
the theoretical predictions, without any parameter tun- 
ing. This provides very strong support for the hydrody- 
namic model description of the fireball evolution, from 
a new set of observables that is quite independent of all 
previously studied observables (pr-spectra, anisotropic 
flow coefficients v n , and HBT radii). 

We note that the non-linear mode-coupling first discov- 
ered in [Io| , and the event-plane rotations driven by this 
non-linear effect, are key to the qualitative agreement be- 
tween theory and data in Figs. [5] and SJ We doubt that a 
similar agreement can be obtained with dynamical mod- 
els that do not rely on a large degree of local thcrmaliza- 
tion in the expanding fireball, or from an approach based 
on linear [7|, [22|, |35| hydrodynamic response to the initial- 
state density fluctuations. Inclusion of first-order non- 
linear terms in the hydrodynamic response Hi| appears 
to yield event-plane correlations with qualitatively simi- 
lar features as shown here [37| . but quantitative success 
likely requires a numerical approach that fully accounts 



2] agree qualitatively, and even almost quantitatively, with the 
corresponding correlations between the angles associated with 
the harmonic coefficients of the anisotropic hydrodynamic flow 
velocity profile along the freeze-out surface. 
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FIG. 3: (Color online) Similar to Fig. [T] but for selected [27| three-plane correlators of the form (cos(ciZ$i+c„n$ n -|-c m m$ m )), 
where the a are integers satisfying cil+c n n+c m m—0 [2^ ]. 
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FIG. 4: (Color online) The final-state event-plane correlators corresponding to the initial-state correlators between three 
participant planes of different harmonic order shown in Fig. [3] Solid and dashed lines show results from viscous hydrodynamics 
with MC-Glauber and MC-KLN initial conditions, evolved with ?7/s=0.08 and 0.2, respectively; filled circles show ATLAS data 
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for the intrinsic nonlincarity of viscous hydrodynamics. 

A closer look at Figs. [1] and [3] shows that the MC-KLN 
model tends to produce stronger correlations between the 
initial-state participant-plane angles $„ than the MC- 
Glauber model. We observe that hydrodynamic evolu- 
tion translates the stronger initial-state participant-angle 
correlations into stronger final-state event-plane correla- 



tions, even though the signs of some of the correlators 
featuring the strongest correlation strengths flip between 
initial and final state. This is especially true for the two- 
plane correlations shown in Fig. [TJ while the three-plane 
correlators exhibit some exceptions to this "rule" in the 
most peripheral collisions. The experimental data ap- 
pear to prefer the stronger angle correlations in the ini- 
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tial profiles from the MC-KLN model, even though this 
model gives an elliptic-to-triangular flow ratio V2/V3 that 
is much larger than measured (lEj , caused by a larger 
£2/23 ratio than in the MC-Glauber model [l0|. These 
observations show that a combined analysis of both the 
anisotropic flow coefficients v n and their associated flow 
angles (and the correlations among them) promises 
to yield powerful constraints on initial state models for 
the fireball energy density profiles created in heavy-ion 
collisions. 

Even though more detailed studies will be necessary to 
fully explore the event-plane correlations discussed in this 
paper, the calculations presented here suggest that very 
likely neither the MC-Glauber nor the MC-KLN initial 
conditions will ultimately provide a quantitatively satis- 



factory description of the experimental data from the AT- 
LAS Collaboration [27j . How to turn the multitude of al- 
ready measured and in the future measurable anisotropic 
flow observables (magnitudes and angles) into a focused 
search for the correct initial-state model is an interesting 
and welcome new challenge for the theory community. 
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